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In this paper, we examine the effects of radiation pressure, 
Poynting-Robertson (PR) drag, and solar wind drag on dust grains 
trapped in mean motion resonances with the Sun and Jupiter 
in the restricted (negligible dust mass) three-body problem. We 
especially examine the evolution of dust grains in the 1 : 1 reso- 
nance. As a first step, the Sun and Jupiter are idealized to both be 
in circular orbit about a common center of mass (circular restricted 
three-body problem). From the equation of motion of the dust 
particle in the rotating reference frame, the drag-induced time rate 
of change of its Jacobi “constant, ” C, is then derived and expressed 
in spherical coordinates. This new mathematical expression in 
spherical coordinates shows that C, in the 1 : 1 resonance, both 
oscillates and secularly increases with increasing time. The new 
expression gives rise to an easy understanding of how an orbit 
evolves when the radiation force and solar wind drag are included. 
All dust grain orbits are unstable in time when PR and solar wind 
drag are included in the Sun-Jupiter-dust system. Tadpole orbits 
evolve into horseshoe orbits; and these orbits continuously expand 
in size to lead to close encounters with Jupiter. Permanent trapping 
is impossible. Orbital evolutions of a dust grain trapped in the 1 : 1 
resonance in the planar circular, an inclined case, an eccentric 
case, and the actual Sun-Jupiter case are numerically simulated 
and compared with each other and show grossly similar time behav- 
ior. Resonances other than 1 : 1 are also explored with the new 
expression. Stable exterior resonance trapping may be possible 
under certain conditions. One necessary condition for such a trap 
is derived. Trapping in interior resonances is shown to be always 

unstable. O 1995 Academic Prea, Inc. 


1. INTRODUCTION 

In a previous paper (Liou and Zook, 1995, hereinafter 
referred to as Paper I) we described the existence and 
evolution of a dust ring, of main belt asteroidal origin, 


composed of micrometer-sized particles that were trapped 
into the 1 : 1 mean motion resonance with Jupiter. It was 
found that a limited size range (radius ~ 1 /xm) of asteroidal 
dust particles, due to the radiation pressure force felt 
when particles are released from their parent bodies, often 
would be injected directly from main belt parent asteroids 
into the 1 : 1 mean motion resonance zone with Jupiter. 
There they would remain trapped for tens of thousands of 
years before having close encounters with Jupiter which 
removed them from the resonance. 

When a particle is trapped, it usually starts with a tad- 
pole-type orbit oscillating around L4 or L5. The maximum 
range of its libration angle increases with time. This angle 
gradually increases to such an extent that the particle’s 
orbit becomes a horseshoe. The horseshoe orbit keeps 
expanding until it becomes possible for the particle to 
have a close encounter with Jupiter and, eventually, the 
particle is scattered out of the resonance. Radiation pres- 
sure, Poynting-Robertson (PR) drag, and solar wind drag 
control the orbital evolution of such particles. Otherwise, 
they would remain trapped in a stable state and behave 
just like the Trojan asteroids. 

The effect of radiation pressure and PR drag at the 
1 : 1 mean motion resonance in the restricted three-body 
problem has been studied to some extent previously. 
These studies include, for example, those of Colombo et 
al. (1966), Chemikov (1970), and Schuerman (1980) who 
discuss the positions as well as the stability of the Lagran- 
gian equilibrium points when radiation pressure and PR 
drag forces are included. More recently, Murray (1994) 
systematically discusses the dynamical effect of general 
drag (nebular drag, gas drag, and PR drag) in the planar, 
circular, restricted three-body problem, including also a 
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discussion of positions and stability of the Lagrangian 
equilibrium points. However, these stability analyses are 
all based on a linear approximation that assumes the dis- 
placements of orbits from the Lagrangian equilibrium 
points are small and that the drag force is also small. 
Secure conclusions could not, from these analyses, easily 
be drawn about the global stability of orbits in the 1:1 
resonance with large displacements from the Lagrangian 
equilibrium points (see, for example, Murray 1994. p. 
480). 

Earlier analyses of the time variation of the Jacobi 
“constant” were used, for example, to study the “Jacobi 
capture” of dust grains near the Earth (Colombo et al. 
1966), to determine the stability of an orbit with a drag 
force proportional to the velocity of a particle in the rotat- 
ing reference frame (e.g., Murray 1994), and the stability 
of Trojan precursors with nebular gas drag (Peale 1993). 
In Murray’s paper (1994), the method is also applied, in 
a two-dimensional analysis, to a simple drag force propor- 
tional to the velocity of a particle in the inertial reference 
frame. He was not able, however, to determine (in Carte- 
sian coordinates) the ultimate stability of the orbit at large 
displacements from the Lagrangian points. 

Building on previous work (especially that of Murray 
1994), we derive a new, and simple, mathematical expres- 
sion in spherical coordinates that accurately describes 
how radiation pressure, PR, and solar wind drag forces 
change the Jacobi constant with time and thereby control 
the orbital evolution of dust particles in a 1 : 1 mean motion 
resonance location. The formula works equally well for 
large displacements as for small. We also apply this new 
expression to the orbital evolution of dust grains trapped 
in other types of resonances. 

With this new mathematical expression, we are also 
better equipped than before to understand the physics 
behind the global evolution of particles, how their orbits 
expand, how their orbits change from tadpoles to horse- 
shoes, etc. We show that 1 : 1 resonance traps are always 
unstable under these drag forces. Tadpole orbits will 
evolve into horseshoe orbits and, eventually, close en- 
counters with Jupiter become possible and the particles 
exit the resonance trap. Hence, no permanent trapping 
is possible for dust grains in the 1 : 1 mean motion reso- 
nance. We also show that interior resonance traps are 
similarly unstable, but that quasi-stable exterior reso- 
nance traps may exist. One necessary condition for quasi- 
stable exterior resonance traps is derived. 

The nature of the radiation force and solar wind drag 
is briefly summarized in Section 2. In Section 3, we derive 
the time dependence of the Jacobi constant on radiation 
pressure, PR drag, and solar wind drag in the circular, 
restricted three-body problem. We also show how the 
time rate of change of the Jacobi constant relates to the 
orbital evolution of particles. The case where all three 


bodies are in the same plane is discussed in additional 
detail. In Section 4, several numerical experiments, in- 
cluding those of noncircular, nonplanar cases are pre- 
sented. The results are compared with the theory from 
Section 3. The implications for micrometer-sized asteroi- 
dal dust particles and Trojan dust particles trapped in the 
1 : 1 mean motion resonance with Jupiter is also discussed. 

2. RADIATION FORCE AND SOLAR WIND DRAG 

The orbital evolution of dust grains in the Solar System 
is affected by the radial solar radiation pressure force, by 
gravitational forces due to the Sun and planets, and by 
PR and solar wind drag forces. In an inertial reference 
frame, the equation of motion of a dust particle with geo- 
metric cross section A and mass m moving under the 
influence of the gravitational forces of the Sun, nine plan- 
ets, the radiation pressure force, and the drag forces can 
be written as (e.g., Burns et al. 1979) 


V = - r , - i + (SA/c)Q pr 

L r s J n = 1 r n 

x[(,_ ( , + sw , l-E.),,- 


( 1 ) 


(1 + sw)- 
c 


where v is the velocity vector of the particle, r s is its unit 
position vector with respect to the Sun, and M 0 is the 
Sun’s mass, while M n and r n are the mass of the nth planet 
and the position vector of the particle with respect to that 
planet. 5, c, and Q pr are the solar energy flux density, 
the speed of light, and the radiation pressure coefficient, 
respectively. The velocity-independent part of the third 
term on the right-hand side of Eq. (1) is the radiation 
pressure force term, while the velocity-dependent parts 
are the drag terms. Solar wind drag is caused by the 
interaction between solar wind ions and the dust particle, 
and the net effect is similar to that of PR drag. The ratio 
of solar wind drag to PR drag is represented by sw in Eq. 
(1). We have assumed sw to be 0.35 (Gustafson 1994). 
The radiation pressure coefficient, Q pr , depends on the 
properties (density, shape, size, etc.) of the particle. Tra- 
ditionally, a dimensionless quantity, /3, is introduced to 
specify the effect of radiation pressure and PR drag. It is 
defined as 


_ radiation pressure force _ SAQ pr r) 
p solar gravitation force cGM Q m * 

The radiation pressure force is the reason some asteroi- 
dal dust particles, with /3 ~ 0.26, are injected into the 1 : 1 
mean motion resonance with Jupiter (see Paper I); it is 
responsible for changing the effective mass of the Sun as 
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seen by the dust particles. Consequently, it changes the 
new location of the 1 : 1 mean motion resonance to a n - 
(1 - P) m a j9 where a } is the semimajor axis of Jupiter’s 
(or any arbitrary planet’s) orbit. PR and solar wind drag 
forces, on the other hand, cause dust particles to lose 
orbital energy and angular momentum and spiral toward 
the Sun. A detailed description of the radiation pressure 
force, PR, and solar wind drag forces can be found in 
Burns et al. (1979). 

3. THE RADIAL RADIATION PRESSURE FORCE, PR, 
AND SOLAR WIND DRAG IN THE RESTRICTED THREE- 
BODY PROBLEM 

3.1. Time Variation of the Jacobi Constant 

In the classical circular restricted three-body problem, 
the three colinear Lagrangian equilibrium points, LI, L2, 
and L3, are unstable while the two equilateral Lagrangian 
equilibrium points, L4 and L5, are stable (with the condi- 
tion that the mass ratio of the secondary to the primary 
is less than 0.0385); this is despite the fact that the L4 
and L5 points are at effective potential maxima in the (x, 
y) plane in the rotating reference frame (e.g., Greenberg 
and Davis 1978). The existence of a dissipative medium 
may destroy the stability of the L4 and L5 regions, how- 
ever, depending on the nature of the dissipative mecha- 
nism (Greenberg 1978, Yoder 1979). For dust particles in 
the Solar System, PR and the solar wind drag serve to take 
away orbital energy and angular momentum expressed in 
an inertial system. Their dissipative nature is well defined 
by the velocity-dependent terms in Eq. (1). These dissipa- 
tive forces on dust grains can sometimes be counterbal- 
anced, at least in part, by planetary gravitational perturba- 
tions, especially when the dust grain is in orbital 
resonance with a planet. As already achieved by others 
(e.g. , Murray 1994), and for the convenience of the reader, 
we first derive the three-dimensional equations of motion 
of the particle in Cartesian coordinates — both in an iner- 
tial reference frame and in a reference frame that rotates 
with the two massive bodies. We then derive the dissipa- 
tive mechanisms of the drag forces in spherical coordi- 
nates in the rotating reference frame and study their ef- 
fects on the stability of the equilibrium points and on the 
overall evolution of an orbit. 

We first define the symbols and units used in the circu- 
lar, restricted three-body problem. Let S and J be repre- 
sented by two masses fi { and /x 2 , respectively (where 
(jl x + M2 = 1» Mi > M 2)* moving in circular orbit around 
their center of mass (CM) due to mutual gravitational 
attraction. The inertial reference frame, (£, 17, £), has its 
origin at CM. (f , 17) is the orbital plane of S and J. The 
rotating reference frame (jc, y, z) is rotating with constant 
angular velocity around the collocated £ and z axes, with 
the x axis lying along the line from CM to J. The distances 


V 



FIG. 1. The inertial coordinate system, (£, tj), and the coordinate 
system, (x, y), corotating with S and J. Both systems have their origins 
at the center of mass (CM) of S and J. The polar coordinate system, 
(r, <f>) is also shown. It is related to (. x , y) by x = r cos <f> and y = r sin 
<t>. See text for the definitions of other quantities. 


from the test particle, p, to CM, S, and J, are r, r,, and 
r 2 , respectively. The unit of length is the distance between 
S and J. The unit of mass is the total of S and J. The unit 
of time is the orbital period of J divided by hr. In such 
normalized units, the gravitational constant and the angu- 
lar velocities of S and J with respect to CM are all unity. 
Figure 1 shows the relationships between these quantities 
in the planar case (where £ and z components are zero). 
The components of Eq. (1), with /3 from Eq. (2), are 

i = “(I - 1 ~ M2 + 0 + SW)F pr , { 

r i '2 

v = -(1 - P)v-\ ^—r 1 ~ M 2 + (1 + sw)f pr (3) 
r \ r 2 

£ = -(1 - /3)/*i^— r 1 - + 0 + sw)f pr-£ . 

r l r 2 

These are the components of the equation of motion of 
a dust grain in an inertial reference frame with solar radia- 
tion pressure and PR and solar wind drag terms. Here the 
radiation pressure force has been combined into the first 
terms on the right-hand side of Eq. (3) and the drag compo- 
nents are 

^PR.f = 

g Mi f [( j ~ + to - m>h + « - CiKKj - €0 , U 

" cr i 1 r\ 
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F pr.t) — 


F PR,{ = 


Mi [[(£ - fi)f + (y- V\)y + (£ ~ i\)tKy ~ V\) , • 
2 ~~ + V 


-fi 


Mlf 

cr}\ 


[(£ - f »)f + ft - vi)v + «- £»)£](£ - W 


n 


+ 4 


Using the coordinate transformation, 

— jc cos t — y sin t 
7) = xsint + y cos r. 


Equation (3) becomes, in the rotating reference frame, 


- (l-0>Mi, , , 

x - 2y = x 3 (x + M2) 


or 


f (X Ml) + (1 + SW)Fp Rjr 
r 2 

y + 2x = y- — — r^-y ~ tjjy + (1 + SW )F PR y 
r, r 2 

2= _<LZ«£! z -^ + ( , + sw) f bu 


All 

x-2y = a -^ + (\ + sw)F PRiJt 

y + 2x = ^ + (1 + sw)F pr 
dy 

z = |^ + 0 + sw)F PR i , 


where 


[/ = 


<iZ»^ + ^ + « (x! + j2) 

r, r, 2 


PR,* 


-0 


Mil 

cr i 1 


[(* + /x 2 )j + (y - /u. 2 )y + zz](x + 


+ (x - y) 


\ [(* + Mz)* + (y- v-i)y + zz]y 


n 


mi f It* + mz)* + (y - Mz)y + zzlz 

1 r? 


^PR.z “ P 2 

Lt 1 


+ (>’ + *) 

+ z(. (9) 


(4) 


When there are no dissipative forces (i.e. , F PR i - 0), there 
is a constant of the motion (e .g. , Moulton 1 9 1 4) , the Jacobi 
constant, defined as 


C ss 2U - (i 2 4- y 2 + z 2 ). 


( 10 ) 


(5) 


If we multiply the first of Eq. (6) by x, the second of Eq. 
(6) by y, and the third of Eq. (6) by z, and add, and 
use Eq. (10), we obtain the time variation of the Jacobi 
constant due to drag forces, 


C = 2(1 + sw)^y (x 2 + y 2 + z 2 )-(yx~xy) 
erf [ 

{xx 4- yy + zz + M 2 x)(xr + yy + zz + fi 2 x - i x 2 y)l 

'i J 


(ID 


( 6 ) 


or 


C = 2(1 + sw)^y (x 2 + y 2 + z 2 ) - (yx - xy) 
erf L 

(xx + yy + zz) 2 


( 12 ) 


+ ■ 


+ C(M2) 


(7) 


Here C(/x 2 ) are terms which are multiplied by the second- 
ary mass, /x 2 , and in the case of the Sun-Jupiter system 
where fi 2 ~ 0.001, these terms can be neglected. This 
formula is similar, in part (the major difference being the 
rf 2 dependence of PR and solar wind drag that we use), 
to that derived in two dimensions by Murray (1994). It 
will prove to be much easier, as we shall see, to visualize 
particle dynamics, and the behavior of C, in a spherical 
coordinate system (r, 0, </>) given by the transformation. 


( 8 ) 


is the modified potential in the rotating reference frame. 
The PR drag components in (x, y , z) coordinate system are 


x = r cos 0 cos </> 
y = rcos 0 sin $ 
Z = r sin 0. 


(13) 


Here r is the radial distance from CM to the particle, 0 
is the latitude angle measured from (x, y) plane toward 
positive z direction, and <f> is the longitude angle measured 
counterclockwise from positive x direction toward posi- 
tive y direction. With this transformation, Eq. (12) be- 
comes 
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C = 2(1 4- sw) 

+ r 2 <b cos 2 0 + 


r 2 + ( rd ) 2 + ( r<f> cos 6) 2 

cr i L 

0 ] 


(14) 


If we now use the law of cosines, (r/r,) 2 can be expressed 
as 

fey— 


where i// is the angle between the positive x axis and the 
position vector of the particle. This angle becomes <f> 
when 0 equals to zero. In the Sun-Jupiter-dust system, 
fi 2 ~~ 0.001, r, ~~ r, and Eq. (14) can be written as 


C = 2(1 + sw)^ 
c 



+ (<£ cos 0) 2 + cos 2 0 




(16) 


This equation shows that when PR and solar wind drag 
are included in the circular Sun-Jupiter-dust system, C 
has three positive-definite terms and one term depending 
on <t> (the last term), which depends on the location of the 
resonance. In a 1 : 1 resonance, the sign of <j> oscillates in 
time and its magnitude is less than unity. In an interior 
resonance f is positive-definite while in an exterior reso- 
nance, <t> is negative-definite. As C (see Eq. (10)) defines 
the locus of coordinates where velocity is zero, an increas- 
ing or decreasing C implies increasing or decreasing exclu- 
sion zones. This expression is independent of how far an 
orbit is from any given Lagrangian equilibrium point. It 
can be used to find the stability of an orbit in a particular 
type of resonance. Detailed discussions of this expression 
in various types of resonances are in Sections 3.3 and 3.5. 

In the classical circular restricted three-body problem, 
the C values associated with five Lagrangian equilibrium 
points characterize certain stable orbits that a test particle 
could have. Except for C values at the L4 and L5 points, 
the C values define the “forbidden regions” (where the 
velocity of the test particle would be imaginary). Forbid- 
den regions and regions where motion is possible are 
separated by zero velocity surfaces or, in the planar case, 
by zero velocity curves (ZVCs). For example, a test parti- 
cle with a near circular orbit with a C value between Qj 
and Cl 4 has a tadpole- type orbit. Such a particle can never 
enter a certain region around L4 (or L5, depending on 
where the particle is placed initially). Here C LN is the 
Jacobi constant associated with the zero velocity surface 
that contains the Lagrangian point L N. A test particle 


with a near circular orbit with a C value between C L2 
and C L3 has a horseshoe-type orbit (e.g., Szebehely 1967, 
Danby 1989). 

In the classical circular restricted three-body problem, 
C is a constant. A test particle in a 1 : 1 resonance cannot 
change its orbit from a tadpole to a horseshoe or vice 
versa. This is no longer true when the drag forces are 
included. The radiation pressure force reduces the effec- 
tive mass of the Sun and changes the locations of the 
Lagrangian equilibrium points. The drag forces also have 
some effect on the locations of the Lagrangian equilibrium 
points as indicated by Eqs. (6) and (9). 

It has been shown previously, in the two-dimensional 
case, that, in a 1 : 1 resonance with PR drag, a test particle 
with a small initial displacement from L4 (or L5) will 
oscillate around that equilibrium point, with the displace- 
ment growing exponentially in time (e.g., Colombo et al. 
1966, Schuerman 1980, Murray 1994). However, when the 
displacement becomes large, their linear approximation 
cannot be used and firm conclusions could not easily be 
drawn about the evolution of an orbit. Our new way of 
using the Jacobi constant, however, has no such limita- 
tion, if the variation in C can be mathematically expressed 
as is the case with Eq. (16). Together with Eq. (10), they 
can be used to determine the time dependence of the 
global behavior of an orbit regardless of how far the orbit 
is from any one given Lagrangian equilibrium point. 

In the Solar System, most of the asteroids are concen- 
trated toward the ecliptic. This is also true for the microm- 
eter-sized asteroidal dust particles that become trapped 
in the 1:1 mean motion resonance with Jupiter, where 
their inclinations are less than 20° (Paper I). In such a 
case, a coplanar circular restricted Sun-Jupiter-dust for- 
mulation is a good first approximation and provides a 
good foundation for understanding dust motion in the 
actual Sun-Jupiter system. Likewise for many dust grains 
originating from Trojan asteroids. In such a case, Eq. (16) 
can be simplified to 


C = 2(1 + sw) 




+ <f> 2 + <j> 


(17) 


3.2. Locations of the Lagrangian Equilibrium Points 
and the Zero-Velocity Curves 

In order to calculate the locations of the Lagrangian 
equilibrium points, Eq. (6) must be solved with the condi- 
tion that all time derivatives are zero. The equations 
become 


0 = x 


(1 ~ ft)M i . M2, , 

3 (* ~1~ M 2 ) 3 (* Ml) 

r l r 2 


, (1 + sw)/3/2, 

rf 




/j. 2 y(x + n 2 ) 


cr 


2 

l 



( 18 ) 
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„ r. o _ ^)^i . 0 + s ' 

°”T — ? — T|J + — ? 

„ f(l — /3 )Mi , M 2 

0 = 1 L — R — + Ti 


sw)/3/x, 

2 

1 


(1 + sw)Pn t n 2 y 


There are two possible solutions in Eq. (20), 


(1-/3)mi M2 (1 + sw)pix lf i 2 y 
r] t cr i 


or 


z = 0. 


The first solution corresponds to equilibrium points out- 
side the (x, y ) plane while the second solution corresponds 
to equilibrium points in the (x, y) plane. We examine first 
the out of the (x, y) plane possibility. It is obvious from 
Eq. (21) that if is no larger than 1, the only possible 
solution is to have y > 0. If we substitute Eq. (21) into 
Eq. (19), we get 


0 = 1 


(1 + SW )(3fl ] X 



cr]y 


(23) 


X (AU) 

FIG. 2. The locations of the five Lagrangian equilibrium points for 
different p values (sw is always 0.35 in our calculation). They are pre- 
sented by black dots. A “O’ 1 indicates the position when p = 0 while 
a “9“ indicates the position when /3 = 0.9. The dots in between “0“ 
and “9“ are each different by 0.1 in p value. The positions of the Sun 
and Jupiter are indicated by crosses. Two thick curves indicate the 
paths of L4 (that in the first quadrant) and L5 (that in the fourth quadrant) 
from p = 0 to p = 0.9 if the effect of drag is ignored; this is the same 
curve as shown by Schuerman (1980) as part of an arc of a circle of 
unit radius centered on the secondary mass (Jupiter). 


With the condition that y > 0, the only possible solution 
is to have x > 0. 

If we now divide both sides of Eq. (18) by (* + nj and 
again use Eq. (21), we get 

0-^- + a[-!-l. < l V W> ^: ? . (24) 

x + /x 2 r{ |_X + fi 2 J crf(x + M2) 


Fig. 2. It can be seen that the radiation pressure force 
has the most effect on changing the equilibrium point 
locations. 

The Jacobi constant associated with the zero-velocity 
curve that contains a given Lagrangian equilibrium point 
can be obtained using Eqs. (8) and (10). The results for 

== 0.26 (or radius ~ 1 fim ), sw = 0.35 particles in a 
circular Sun-Jupiter-dust system are listed in Table I. 

All the Jacobi constants in Table I are less than 3. This 


With the conditions that x > 0 and y > 0, the right-hand 
side of Eq. (24) can never equal to 0. Therefore, we can 
conclude that, when the radiation pressure force on a 
particle is smaller than the gravitational attraction from 
the Sun, there are no equilibrium points outside the 
(x, y) plane. This is also easy to show geometrically. 

For the z = 0 case, Eqs. (18) and (19) can be solved 
using Taylor series expansions in r x and r 2 and some ap- 
proximations (e.g., Schuerman 1980, Murray 1994) or, 
more directly, can be solved numerically. The latter 
method is used here. We show the locations of the five 
Lagrangian equilibrium points for different (3 values in 


TABLE I 

Locations of the Lagrangian Equilibrium Points and the Jacobi 
Constants Associated with the ZVCs That Contain those Points 



X (AU) 

Y (AU) 


LI 

4.5825 

-1.2E - 4 

2.4688134 

L2 

5.4560 

0.0 

2.5468315 

L3 

-4.7047 

9.64E - 2 

2.4555497 

L4 

2.0866 

4.2128 

2.4539510 

L5 

2.1561 

-4.1776 

2.4539510 


Note, p = 0.26 and sw = 0.35, in the circular Sun-Jupiter-dust system 
(the semimajor axis of Jupiter is set to 5.2 AU). 
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-5 0 5 

X (AU) 

FIG. 3. Five Lagrangian equilibrium points and three critical ZVCs 
of the Sun-Jupiter-dust system with /3 = 0.26, sw = 0.35. The positions 
of the Sun and Jupiter are indicated by crosses. While LI (between 
Jupiter and Sun) and L2 (outside Jupiter, along +x direction) are along 
the x axis, L3 has a noticeable nonzero y component. The thin, intermedi- 
ate, and thick curves are the ZVCs that contain L3, LI, and L2, respec- 
tively. 


is due to the effect of radiation pressure. The positions 
of five Lagrangian equilibrium points and the ZVCs that 
pass through LI, L2, and L3 points in the (jc, y) plane are 
shown in Fig. 3. Apparently the radiation force does not 
change the general shapes of these ZVCs. This can be 
seen by comparing Fig. 3 with the classical ZVCs (e.g., 
Szebehely 1967). However, there are two major differ- 
ences. The first one, as pointed out by Colombo et al. 
(1966) and more clearly demonstrated by Murray (1994) 
and depending on the nature and strength of the drag 
force, is that the locations of LI, L2, and L3 may have 
nonzeroy components. With PR and solar wind drag, the 
existence of nonzero y components for LI, L2, and L3 
can be verified easily from Eq. (19). However, only the 
y components of L3 and LI (with very large f3) are notice- 
able in Fig. 2. The second difference for our case of /3 = 
0.26, contrary to the classical /3 = 0 case, is that C LI is 
smaller than C L2 . Consequently, as C increases from C L3 , 
the corresponding horseshoe orbit expands and evolves 
into either a prograde orbit inside an inner oval (containing 
S), or a retrograde orbit outside an outer oval (contain 
J) when the C value reaches C u . As C further increases, 
the ZVC eventually divides the plane into three possible 
regions for an orbit (when C reaches C L2 ). They are a 


prograde orbit inside a small oval around S, an orbit inside 
a smaller oval around J, and a retrograde orbit outside 
the S-J system. Here the “prograde” and ‘‘retrograde 1 ’ 
orbits all refer to the rotating reference frame. 

3.3 . Evolution of an Orbit in a 1:1 Resonance and the 

Effect of Inclination and Eccentricity 

For a particle trapped in a 1 : 1 mean motion resonance 
in the planar case, the sign of </> oscillates with time and 
its magnitude is less than unity. This means that C varies 
periodically with </> and secularly increases with time due 
to the </> 2 and (r/r) 2 terms (see Eq. (17)). The <£ 2 is the 
dominant secular term for a low eccentricity orbit. This 
secular increase in C leads to the expansion of the forbid- 
den region where particle motion is imaginary. This gradu- 
ally forces the tadpole or horseshoe orbits to expand. For 
example, if a particle starts in near circular orbit with a 
small displacement (“small 11 is not a necessary condition) 
from L4 in the rotating reference frame, its corresponding 
C value will be slightly larger than C L4 . The forbidden 
region then, which contains L4, is small. The particle 
moves clockwise around, and exterior to the tadpole- 
type forbidden region. As the particle moves around the 
forbidden region, C undergoes a sinusoidal variation due 
to </> and a smaller secular increase due primarily to <j> 2 . 
This increase in C causes the forbidden region to expand. 
This, in turn, forces the particle’s orbit to expand. When 
C reaches C L3 , the forbidden region becomes a horseshoe, 
and the particle must orbit outside that horseshoe. As C 
further increases, the horseshoe orbit expands to such 
an extent that close encounters with J become possible. 
Eventually, a close encounter with J will remove the parti- 
cle from the resonance. Thus, all orbits are unstable. No 
permanent trapping in the 1 : 1 mean motion resonance is 
possible unless a particle is put right at a given Lagrangian 
equilibrium point, which is highly unlikely in the real Solar 
System (Paper I). 

In the nonplanar case, where the orbit of the dust parti- 
cle is not in the orbital plane of the planet, the conclusion 
that all orbits are unstable is still true as shown by Eq. 
(16). In this case, it is still the expansion of the zero- 
velocity surface that forces the orbit to expand. However, 
because of the effect of 0 on both <j> and </> 2 , and the 
additional term of if, the short-term variation as well 
as the secular increase components will be somewhat 
different from those of the planar case. In addition to the 
projection effect of on the (jc, y) plane, the amplitude 
of the short-term variation will be smaller by a factor of 
cos 2 0. The net result depends on both the inclination and 
the orientation of the orbit (i.e., where the longitude of 
the ascending node is). 

In the classical elliptic restricted three-body problem 
where the planet is in an eccentric orbit and where no 
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radiation force is involved, the Jacobi integral is no longer 
a constant of the motion. The Jacobi constant becomes 
a function of the eccentricity of the planet and time (e.g., 
Szebehely 1967, Marchal 1990). The positions of the La- 
grangian equilibrium points also vary. Here we limit our- 
selves to a simple qualitative discussion on the effect of 
eccentricities of dust grain and of Jupiter’s orbits. When 
the eccentricities of both Jupiter and the dust particle are 
included in the system, additional time-varying periodic 
factors will be combined with r, c/>, and 6 in Eq. (16) to 
modify the problem. This does not change the general 
conclusion that C has a large short-term variation and 
small secular increase components. All orbits are unstable 
under PR and solar wind drag. However, because of the 
strong effect of the test particle’s eccentricity on its r, the 
overall evolution will be sped up and the speed up rate 
increases as the test particle’s eccentricity increases. 

All the above analytical results are verified in our nu- 
merical simulations in Section 4. 

3.4. Damping of the Inclination 

PR drag is an in-plane drag force and has no effect on 
the inclination of an orbit in a Sun-dust two-body system 
(e.g., Burns et al. 1979). However, in a Sun-planet-dust 
three-body system, the inclination of the dust grain can 
be changed. This is especially obvious when the dust grain 
is in a 1 : 1 mean motion resonance with a planet. If we 
substitute from Eq. (9) into the z equation in Eq. 
(6), we have 


(1 - j3)/xj fi 2 


-[■ 

-[(l + SW ) ^]i, 


+ — f + (I + sw) 


y8/x,(rr) 

erf J 


(25) 


where the /x 2 terms in z have been dropped. This equa- 
tion describes a typical ‘‘underdamped” oscillator (e.g., 
Marion 1970), with an e-folding time 


toward and become trapped into an interior resonance 
with a planet, although it has been argued analytically, 
using Hamiltonian mechanics, to be impossible (e.g., 
Henrard and Lemaitre 1983, Peale 1986), such trapping 
has actually been observed in numerical simulations 
(Jackson and Zook 1992, Liou 1993, Marzari and Vanzani 
1994). In all numerical simulations, the trapping is found to 
be only temporary. Expression Eq. (16) does not predict 
whether or not trapping into interior resonances is possible. 
However, it does show that trapping intervals in interior 
resonances are always only temporary. 

If a dust particle does get trapped in an interior reso- 
nance with a planet, its orbital angular velocity is larger 
than that of the planet. Thus, in the rotating reference 
frame, the particle’s <£ is positive all the time. This implies 
all terms on the right-hand side of Eq. (16) are positive. 
As C increases, without limit, the dust grain is eventually 
forced to leave the resonance and continue its journey 
toward the Sun. Permanent trapping will only be possible 
if a case can be found where C, on average, is zero. 

In an exterior resonance, the particle’s </>, in the rotating 
reference frame, is negative all the time. Therefore, in 
order to reach a possible equilibrium trapping, ( Hr ) 2 and/ 
or B 2 and (</> cos B) 2 must just be large enough to counter- 
balance the effect of <j> and lead to C = 0. This implies 
that in order to reach an equilibrium exterior trap, the 
eccentricity and/or inclination of the dust grain must 
reach certain values. Planetary resonant perturbations 
(e.g., Weidenschilling and Jackson 1993) then are required 
to keep the orbital eccentricity at the critical value of e 
for C to equal zero, in the planar case. 

Another way of treating this problem is to look at the 
variation ofC in orbital parameter space. In the planar case, 



'•-(TTSte- (26) 

where all quantities are in the normalized units. In a 1 : 1 
mean motion resonance, the damping in the amplitude of 
the z motion leads to the damping in the inclination of the 
dust grain. The e-folding time here is much longer than 
the PR drag life time of a dust grain, if it is not trapped 
in any resonance. 


The second bracket on the right-hand side is zero and the 
rates of change in a and e due to PR drag are (Wyatt and 
Whipple 1950) 


da\ _ -/3GAf 0 (2 + 3e 2 ) 
,dt) pr c a(l - e 2 ) V2 


and 


3.5. Applications to Interior and Exterior Resonances 

Our result (Eq. (16)) can be applied to study interior \—\ - j e 

and exterior resonances as well. For a dust grain to drift Wf /pr c 2a\\ — e 2 ) U2 


( 29 ) 
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The expression of C in terms of the orbital elements of 
an orbit is the Tisserand criterion (e.g., Moulton 1914) 

c = (* ~ ft) + 2V(1 - /3)a(l - g 2 ), (30) 

a 


where a and e are the dust grain’s semimajor axis and 
eccentricity, respectively. This expression is valid when 
the planet is in circular orbit and when the dust grain is 
not too close to the planet. Even though the dust grain 
can come close to the planet during part of its revolution, 
this expression is still good in describing the general be- 
havior of the dust grain’s orbit. If we calculate the partial 
derivatives of C with a , and e y using Eq. (30) and substitute 
the results, together with Eqs. (28) and (29), into Eq. (27), 
the condition for C = 0 is 



2 + 3r> 2 
2(1 -e 2 ) V2 ‘ 


(31) 


When a dust grain is in a (p + q)ip exterior mean motion 
resonance with a planet, 


a 3 ] l/2 + g 

_l“/3j p 


(32) 


where p and q are positive integers (see also Weidenschil- 
ling and Jackson 1993). This leads to 


2 + 3e 2 _p + q 
2(\-e 2 ) V2 ~ p * 


(33) 


This means, for a given (p + q)lp exterior mean motion 
resonance, there exists an eccentricity (of the dust grain) 
according to Eq. (33) that ensures C = 0. It is not yet clear 
whether completely stable external resonances actually 
exist, but C = 0 is a necessary condition if such stable 
resonances exist. We have noted that in long-lived three- 
body resonances, the eccentricity of the dust grain drifts 
close to the value of e obtained from Eq. (33). In fact the 
e obtained from Eq. (33) seems to give a better limiting 
value of e than the e max given by Eq. (18) in Weidenschil- 
ling and Jackson (1993). 

In a three dimensional case, however, the condition for 
a stable trap may not exist, due to the effect of drag on 
the orbital inclination. 


4. NUMERICAL EXPERIMENTS 

4.1 . 1:1 Mean Motion Resonance in the Coplanar 
Circular Sun-Jupiter-Dust Case 

In this section, we show the results from our numerical 
simulations. The equations of motion of the Sun, Jupiter, 
and a dust particle are integrated using RADAU (Everhart 
1985) on a HP9000 715/75 workstation. The Sun and Jupi- 


4.4 



FIG. 4. Positions of the test particle, with = 0.26 and sw = 0.35 
in the rotating reference frame for the first 300 years in the planar, 
circular Sun-Jupiter-dust system. The Lagrangian equilibrium point L4 
is labeled by the black dot. 


ter are gravitationally interacting in circular orbit with 
each other and the dust particle while the dust particle 
does not affect the motions of the Sun and Jupiter. Jupi- 
ter’s semimajor axis is set to 5.2 AU. The eccentricity of 
the test particle is initially set to near zero to examine 
the circular orbit case. We use fi = 0.26 and sw = 0.35, 
corresponding to a l-/xm-radius dust particle of density 
2.7 gm cm -3 trapped in the 1 : 1 mean motion resonance 
with Jupiter (Paper I). 

The first test particle is initially placed very close to 
L4. It is trapped in the 1 : 1 mean motion resonance with 
Jupiter for several hundred thousand years before it has 
a close encounter with Jupiter at 250,651 years and is 
removed from the resonance. Initially the particle is in a 
tadpole orbit around L4 (Fig. 4). Its orbit expands gradu- 
ally and becomes a horseshoe after 234,200 years (Fig. 
5). The variation in C over the particle’s trapping life time 
is shown in Fig. 6. Note that this particle’s orbit becomes 
a horseshoe when its Jacobi constant, C, exceeds C L3 
(Table I), as we expected from the analysis in Section 3. 
Although this causes the oscillatory amplitude of to 
suddenly increase, and therefore C merely continue to 
smoothly increase with time. 

In Fig. 7, we show the detailed variation of C from 
200,000 to 204,000 years, as an example to demonstrate 
Eq. (17). Qualitatively, it agrees perfectly with Eq. (17), 
with a large short-term variation and a small secular in- 
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X (AU) 

FIG. 5. Positions of the same test particle from Fig. 4, in the rotating 
reference frame at four different time periods. They are, from inside 
out, 140,000 to 141,000 years, 200,000 to 201,000 years, 230,000 to 
231,000 years, and 246,000 to 247,000 years. We connect positions that 
are two years apart to form the curves. The five Lagrangian equilibrium 
points are indicated as black dots. The positions of Sun and Jupiter are 
labeled by crosses. 


crease. The time derivative of C from Fig. 7 is also com- 
pared, quantitatively, with the right-hand side of Eq. (17). 
The comparison is in Fig. 8a. The curve is the time deriva- 
tive of C from Fig. 7 while the points are the results by 
using positions of the test particle as functions of time to 
calculate the right-hand side of Eq. (17). The agreement 
is excellent. Note the large oscillatory variation, due to 
</>, and the smaller positive secular increase with time, 
mostly due to <j> 2 (Fig. 8b). Figure 8c shows the still smaller 
secular increase due to (r/r) 2 . These three figures show 
that Eq. (17) does describe the time variation of the Jacobi 
constant correctly. 

If a test particle is placed initially near L5 instead of 
L4, its evolution is similar to the above case. Based on 
our analysis in Section 3, there is no difference between 
C LA and C L5 for particles with f3 as large as 0.26, and 
sw = 0.35, in the restricted Sun-Jupiter-dust system. 
Our numerical experiments also show no difference in 
terms of the global behavior between these two cases. 

42. 1 : 1 Mean Motion Resonance in the Nonplanar 

and Noncircular Sun-Jupiter-Dust Cases 

In this section, we show some numerical examples of 
the nonplanar, noncircular cases. The first one is a dust 
particle that starts with an eccentricity of 0.1 and the 
longitude of perihelion at 0°, moving in the Sun-Jupiter 
system. Jupiter is assumed to have a circular orbit. This 
is still a circular (circular Jupiter orbit) case, but the eccen- 



FIG. 6. The Jacobi constant of the test particle over a period of 250,651 years in the coplanar, circular restricted Sun-Jupiter-dust system. 
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FIG. 7. A detailed plot of part of Fig. 6, from 200,000 to 204,000 years. 


tricity of the dust particle, which decays with time, makes 
the overall evolution quite different from the previous 
case. The dust particle is initially placed near L4. It starts 
with a tadpole orbit and the orbit gradually expands and 
becomes a horseshoe at 193,660 years (Fig. 9). The horse- 
shoe orbit keeps expanding until the particle has a close 
encounter with Jupiter at 204,476 years. The particle is 
then removed from the resonance. In Fig. 10 we show 
the variation of C of the particle. Detailed examination 
of the curve shows that the particle’s orbit becomes a 
horseshoe when its C value exceeds C L3 . The short-term 
variation does not seem to be too much larger than that 
in the planar, circular case, while the secular increase is 
more dramatic, due to the particle’s eccentricity effect 
on r. The total trapping life time is significantly shorter 
than that of the circular case. 

The next example is a dust particle with 5° inclination 
with respect to Jupiter’s orbital plane. The longitude of 
the ascending node is 30°. Both Jupiter and the dust parti- 
cle have circular orbits. This dust particle starts its orbit 
as a tadpole around L4. Its orbit expands and becomes 
a horseshoe at 236,656 years. Finally the particle has a 
close encounter with Jupiter at 245,579 years and is re- 
moved from the resonance. The variation in C of the 
particle is shown in Fig. 11. Qualitatively, its overall be- 
havior is similar to the circular planar case with a large 
short-term variation; however the initial secular increase 
is faster with time because of the 6 2 components. This 


fast increase slows down as the dissipative nature of drag 
contracts the particle into the (. x , y ) plane (i.e., reduces 
its orbital inclination) where the potential is a minimum 
along the z direction. The e-folding time according to Eq. 
(26) is about 202,500 years, which is in good agreement 
with the decay in inclination from the numerical result 
(Fig. 12). There is no significant change in the total trap- 
ping life time for this inclination compared with the planar 
case. 

A third example is a dust particle with eccentricity of 
0.1, longitude of perihelion at 70° and inclination of 5° 
(with respect to the orbital plane of Jupiter) and longitude 
of the ascending node at 30°. The eccentricity of Jupiter 
is 0.048. This particle starts with a tadpole orbit around 
L4 and becomes a horseshoe at 223,600 years. It has a 
close encounter with Jupiter at 229,463 years and is re- 
moved from the resonance after that. The variation in C 
of the particle is shown in Fig. 13. Both its long-term and 
short-term variations are much more complicated than 
the previous cases. However, it is still the secular increase 
component that gradually enlarges the forbidden regions 
that forces the orbit to expand from a tadpole to a horse- 
shoe and to have a close encounter with Jupiter. 

As we can see from the above examples, the evolution 
of an orbit becomes more complicated when the eccentric- 
ity and inclination of the dust grain are involved. Likewise 
when Jupiter’s orbital eccentricity is considered. How- 
ever, although the planar circular analysis in Section 3 
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FIG. 8. (a) Comparison between the time derivative of C from Fig. 

7 (shown as curve) and the results (shown as black dots) by using 
positions of the dust particle as functions of time to calculate the right- 
hand side of Eq. (17). The agreement is excellent, (b) The contribution 
of <f> 2 term in the bracket of Eq. (17) to the secular increase of C. Term 
2 = 2(1 + sw)(/3/i|/c)<^ 2 . (c) The contribution of 2(r/r) 2 term in the bracket 
of Eq. (17) to the secular increase of C. Term 1 = 2(1 + sw)(/3^ ( /c)[2(f/ 
r) 2 ]. 


does not describe in detail the evolution of an orbit in these 
cases (especially the last one), it still provides guidance in 
understanding the nature of the evolution. The natural 
trend for an orbit under the influence of PR and solar 
wind drag is to evolve from a tadpole to a horseshoe 
which eventually leads to a close encounter with Jupiter. 
An analytical method for the eccentric Jupiter system 
with mutual inclination between objects may be possible 
to develop, but much can be learned from direct numerical 
simulation alone. The basic physics involved is the same 
as in the case of the planar circular system. 

4.3. Implication to Asteroidal Dust Particles and 
Trojan Dust Particles Trapped in the l : l Mean 
Motion Resonance with Jupiter 

The trapped ring particles we predicted to exist from our 
numerical analyses in Paper I were of main belt asteroidal 


origin; they had £ - 0.26 and sw = 0.35 and were often 
injected right into the 1 : 1 mean motion resonance loca- 
tion (which, with this value of /3, is somewhat interior 
to Jupiter’s orbit) when they were released from their 
asteroidal parent bodies. In addition to Jupiter, there were 
six other planets perturbing their orbits (in our numerical 
simulation, we did not include Mercury and Pluto). The 
initial eccentricities of the trapped dust grains lay between 
0.1 and 0.35 while their typical inclinations lay between 
0° and 20°. Their typical initial libration angles (angle from 
Jupiter) were about 45° to 135 °. While they were trapped, 
their orbits expanded from tadpole orbits to horseshoe 
orbits. Most of the particles got out of the resonance 
within 30,000 years. Close encounters with Jupiter caused 
them to be removed from the resonance (Paper I). Dust 
particles, with considerably smaller /3 values, produced 
from Trojan asteroids can also be trapped in the 1 : 1 
resonance (Paper I). Their orbital evolution is similar to 
that of the asteroidal dust particles. As a matter of fact, 
any dust particles trapped in a 1 : 1 resonance with any 
planet will behave similarly. Their evolution can all be 
described by the analyses in Section 3. 

The real Solar System is more complicated than a sim- 
ple circular restricted Sun-Jupiter-dust system. But the 
physical nature of the evolution of those dust particles is 



X (AU) 

FIG. 9. Positions of a test particle that started with 0.1 eccentricity 
and 0° longitude of perihelion in the rotating reference frame at two 
different time periods. The Sun and Jupiter are in circular orbits and 
the dust particle orbits in the same plane. The inner curve is from 140,000 
to 141,000 years (the dust particle’s eccentricity has dropped to less 
than 0.01). The outer curve is from 200,000 to 201,000 years. Again, 
we connect points that are two years apart to form the curves. 
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FIG. 12. The decay of orbital inclination of the test particle in Fig. 11. The predicted ^-folding time from Eq. (26) is about 202,500 years. 



Time (years) 


FIG. 13. The Jacobi constant of a test particle over the period of 229,463 years. The test particle initially has an inclination of 5° (with respect 
to Jupiter’s orbital plane) with longitude of the ascending node at 30° and an eccentricity of 0.1 with longitude of perihelion at 70°. The eccentricity 
of Jupiter is 0.048. Its longitude of perihelion is set to 0°. 
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still qualitatively described by our analysis. When drag 
forces are negligible, particles can get trapped for a much 
longer time scale, as in the case of the present Trojan 
asteroids. When PR and solar wind drag are included, the 
Jacobi constant of a dust particle will increase as is shown 
by Eq. (16). This leads to the expansion of an orbit from a 
tadpole to a horseshoe and, eventually, a close encounter 
with Jupiter becomes unavoidable. This is in good agree- 
ment with the results from Paper 1. 

We can also compare the trapping life time of dust 
particles in Paper I with our analysis in this paper. In 
the planar, circular Sun-Jupiter-dust system, if a dust 
particle has — 0.26, sw = 0.35 and starts with a libration 
angle from 45° to 135° (similar to the second inner most 
orbit in Fig. 5), it will remain trapped for approximately 
50,000 years before having a close encounter with Jupiter 
(use the data from Section 4.1). If we consider the effects 
of eccentricity (0.048 for Jupiter and up to 0.35 for dust 
particles) and inclination (up to 20°) on the system, they 
have the ability to speed up the growth of an orbit through 
increased 0 2 and ( rlr ) 2 in Eq. (16); with the additional 
perturbation from other planets, 30,000 years seems to 
be a reasonable trapping life time for real asteroidal dust 
particles trapped in the 1 : 1 mean motion resonance with 
Jupiter. 

Another phenomenon about the 1 : 1 trap, from our nu- 
merical calculations in Paper I, is that it is impossible for 
a dust grain to drift toward and get trapped in the 1:1 
resonance with Jupiter. This can easily be understood 
from the analyses in Section 3. When a dust grain is 
approaching, under PR and solar wind drag, from outside 
the location of 1 : 1 resonance toward it, the C of the dust 
grain decreases. When C reaches C L , , trapping in the 1 : 1 
resonance becomes possible. However, if the dust grain 
does get trapped, its C immediately starts to increase (Eq. 
(16)) and expels the dust particle from the 1 : 1 resonance. 
Thus, it is impossible for a dust grain to drift toward and 
get trapped in a 1 : 1 resonance. Actually, when the planet 
involved is a giant planet, such as Jupiter, gravitational 
scattering usually prevents the dust grain from even reach- 
ing C L1 . 


5. CONCLUDING REMARKS 


Instead of the Jacobi constant, C, one can consider the 
“modified total energy “ in the rotating reference frame 
defined by 


£,--§ = |v’-(A 


(34) 


where 


r i (1 -/3)^l M2 1 , , , „ 

-u= — -r(x 2 + /) 

' 1 • 1 £ 


(35) 


is the Jacobi potential. Without drag, E } is also a 
constant of the motion with a modified potential, — (/, 
due to the rotation of the coordinate system (see, e.g., 
Danby 1989). It is also equal to the Hamiltonian constant 
in the rotating system (Szebehely 1967). This constant 
comes from the only explicitly determined integral of 
motion in the three-body system. Because of its close 
resemblance to the energy in the inertial system, it may 
be more appropriate to call E } the Jacobi energy. 

In the two-body Sun-dust problem the effect of the 
drag forces is to take away orbital energy and angular 
momentum from the particle (Wyatt and Whipple 1950, 
Burns et ai 1979). In the three-body rotating reference 
frame, drag acts to take away the modified energy, or 
Jacobi energy, of the particle while the grain is in a 1 : 1 
trap. By following the variation of the Jacobi energy with 
time, we can gain some understanding of the physical 
nature of how an orbit evolves. 

If there were no drag forces, E } like C would remain 
constant with time. At those points where E } = -{/, the 
velocity is zero and we can think of the particle as pressed 
up against the side of the potential energy “hill.” Coriolis 
forces (e.g., Greenberg and Davis 1978) then constrain 
the particle motion to librate around the hill — which may 
be either tadpole-shaped or horseshoe-shaped — leading 
to tadpole- or horseshoe-shaped orbits. 

The PR and solar wind drag forces cause E } to decrease 
with time (equivalent to the increase of C). This, in turn, 
forces the particle farther down the side of the potential 
hill and increases the region of exclusion (i.e., the total 
energy, E } , must always be greater than or equal to the 
potential energy, — U ). Thus PR and solar wind drag 
forces cause the corresponding tadpole or horseshoe or- 
bits to expand with time because dust particles must orbit 
at lower energies about an ever increasing circumference 
of the inaccessible part of the potential hill. 

Relating the circular restricted case to the actual case 
of the Sun-Jupiter-dust system, where Jupiter is in an 
elliptical orbit about the Sun and gravitational perturba- 
tions from other planets must be considered, is still not 
straightforward, however, and we must rely on numerical 
modeling. But we believe the insights that we have gained 
with our new mathematical formulation helps in gaining 
a physical intuition for the cause of the ever increasing 
libration amplitudes. 
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Note added in proof. Since acceptance of our present paper for 
publication, we have become aware of a very interesting paper by 
Beaug* and Ferraz-MeLIo ( Icarus 110, 239-260, 1994). They use a Hamil- 
tonian treatment of the variation of resonance variables to derive their 
Eq. (15), and corresponding Table I, which give relationships between 
orbital eccentricity and integer resonance ratios identical to the require- 
ment we derive in our Eq. (33) for potentially stable trapping in exterior 
resonances. Our treatment differs from theirs in that it depends on an 
examination of the time variation of the Jacobi “constant.” 

REFERENCES 

Burns, J. A., P. L. Lamy, and S. Soter, 1979. Radiation forces on 
small particles in the solar system. Icarus 40, 1-48. 

Chernikov, Y. A. 1970. The photogravitational restricted three-body 
problem. Sov. Astron . 14, 176-181. 

Colombo, G., D. A. Lautman, and 1. 1. Shapiro, 1966. The Earth’s 
dust belt: Fact or fiction? 2. Gravitational focusing and Jacobi capture. 
J. Geophys. Res. 71, 5705-5717. 

Danby, J. M. A. 1989. Fundamentals of Celestial Mechanics , second 
ed., Willmann-Bell, Richmond, VA. 

Everhart, E. 1985. An efficient integrator that uses Gauss-Radau 
spacings. In Dynamics of Comets: Their Origin and Evolution (A. 
Carusi and G. B. Valsecchi, Eds.), pp. 185-202. Reidel, Boston. 
Greenberg, R. 1978. Orbital resonance in a dissipative medium. Icarus 
33, 62-73. 

Greenberg, R., and D. R. Davis, 1978. Stability at potential maxima: 
The L4 and L5 points of the restricted three-body problem. Am. J. 
Phys. 46, 1068-1070. 

Gustafson, B. A. S. 1994. Physics of zodiacal cloud. Ann. Rev . Earth 
Planet. Sci. 22, 553-595. 

Henrard, J., and A. Lemaitre, 1983. A second fundamental model 
for resonance. Celest. Mech. 30, 197-218. 


Jackson, A. A., and H. A. Zook, 1992. Orbital evolution of dust 
particles from comets and asteroids. Icarus 97, 70-84. 

Liou, J. C. 1993. Dynamical Evolution of Asteroidal and Cometary 
Particles and Their Contribution to the Zodiacal Cloud. Ph.D. thesis, 
University of Florida. 

Liou, J. C., and H. A. Zook, 1995. An asteroidal dust ring of micron- 
sized particles trapped in the 1 : 1 mean motion resonance with Jupiter . 
Icarus 113, 403-414. 

Marchal, C. 1990. The Three-body Problem. Elsevier, Amsterdam. 

Marion, J. B. 1970. Classical Dynamics of Particles and Systems , 
second ed. Academic Press, New York. 

Marzari, F., and V. Vanzani, V. 1994. Dynamical evolution of inter- 
planetary dust particles. Astron. Astrophys . 283, 275-286. 

Moulton, F. R. 1914. An Introduction to Celestial Mechanics , second 
ed. Dover, New York. 

Murray, C. D. 1994. The dynamical effects of drag in the circular 
restricted three-body problem: I. The location and stability of the 
Lagrangian equilibrium points. Icarus 112, 465-484. 

Peale, S. J. 1986. Orbital resonances, unusual configurations and exotic 
rotation states among planetary satellites. In Satellites (J. A. Bums 
and M. Matthews, Eds.), pp. 159-223. Univ. of Arizona Press, 
Tucson. 

Peale, S. J. 1993. The effect of the nebula on the Trojan precursors. 
Icarus 106, 308-322. 

Schuerman, D. W. 1980. The restricted three-body problem including 
radiation pressure. Astrophys. J. 238, 337-342. 

Szebehely, V. 1967. Theory of Orbits. Academic Press, New York. 

Weidenschilling, S. J., AND A. A. Jackson, 1993. Orbital resonances 
and Poynting-Robertson drag. Icarus 104, 244-254. 

Wyatt, S. P., and F. L. Whipple, 1950. The Poynting-Robertson 
effect on meteor orbits. Astrophys. J. Ill, 134-141. 

Yoder, C. F. 1979. Notes on the origin of the Trojan asteroids. Icarus 
40, 341-344. 




